Collective effects of the Coulomb interaction in anharmonic quantum dots 
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Deviations from the uniform oscillator spacing, related to the shape of the confining potential, 
have a strong influence on few-electron states in quantum dots when Coulomb effects are included. 
Distinct signatures are found for level spacings increasing (pot shape) and decreasing (cone shape) 
with energy. Cone deformations affect the levels near the ground state, such that observable effects 
are predicted. Pot deformations partially negate the effect of the Coulomb force, thus their spectra 
are similar to those of a perfect oscillator with a smaller Coulomb repulsion. The Coulomb force is 
treated by exact diagonalization, for which purpose efficient closed-form expressions for the matrix 
elements are derived. The Coulomb matrix element in relative coordinates is reduced to a single 
i-O , sum over four binomial coefficients, for which simple analytic approximations are found. 
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I. INTRODUCTION 

Quantum dots are both the smallest electronic devices, and the largest atoms. Technologically, they allow the 
storage and manipulation of single electrons. Physically, they are a practical laboratory to observe few-body Coulomb 
states at energy scales far below those of real light atoms. In small quantum dots, the Coulomb effects are typically 
about one-half of the oscillator spacing, so they have been treated perturbatively with success^ even though they 
would appear a priori to belong to a crossover regime. 

More recent developments have stressed the role of electronic correlations in few-electron dots, among others 
with a view to manipulate quantum states for purposes of computation. Such efforts fall into two categories. One 
aims to construct "quantum dot molecules" in particular geometric patterns^ The other has focussed on "electronic 
''O 1 molecules," referring to the limit of strongly interacting electrons within a single quantum dot, where a molecule-like 
partition of energy into rotational and vibrational degrees of freedom is observed^ These latter states might better 
be called "quantum dot nuclei," since collectivization by the repulsive Coulomb force occurs in the absence of centers 
of opposite charge, which exist in molecules. 

The present paper investigates the particular effect of well shape (in vertical, not horizontal profile), modelled as 
departure from the uniform single-particle energy spacing, on the collectivity of few-electron quantum states within 
a single dot. The modelling is implemented simply by inserting diagonal energy values, other than the "correct" 
oscillator ones. The two qualitatively different cases are level spacings increasing and decreasing with energy, referred 
to respectively as pot-shaped and cone-shaped dots. 

The main result is that cone anharmonicities make quantum dots strongly prone to creation of collective states. The 
physical reason is first demonstrated on the classic case of two electrons ( "quantum dot helium" ) , at modest values 



of the Coulomb force, consistent with the literature. With four electrons, there appear both a change in ordering of 
the lowest excited states, and "intruder" states of exceptionally large spin and orbital momentum among the higher 
t-H . ones. Both of these should lead to observable effects, if such dots can be manufactured. 

Two formal developments are also presented, which are not strictly necessary for the proof of principle given here, 
.£h ! but may be useful in larger calculations. First, an advanced combinatorial technique^ is used to bring the center-of- 
mass (CM) Coulomb matrix elements to a more tractable form. It gives rise to two simple analytical approximations, 
one of them already known^ which give good estimates of the Coulomb matrix elements over a large range of quantum 
numbers. Second, a simple expression for the Coulomb matrix elements in laboratory coordinates is obtained, based 
on an old symmetry-based insight^ The two together mean that binomial coefficients suffice in Coulomb calculations 
for the two-dimensional oscillator: no gamma or hypergeometric functions need apply. 



II. CENTER-OF-MASS MATRIX ELEMENTS 



The 2D oscillator Hamiltonian is 
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where the dimensionless coordinate r has been scaled in units of the oscillator characteristic length Iq = \fhjjj 
It has the eigenfunctions ^ n i, 



vMr) = * nl (x, y) = J , U [ , v L' n (x 2 + y 2 )(x + z^e"^)/ 2 , (2) 
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where L l n is the generalized Laguerre polynomial, and there are 2n + I + 1 degenerate eigenfunctions for each energy 
level: 

Hosc^ni = (2n + l + l)hoj^ n i. (3) 
One can exploit a symmetry of the Laguerre polynomials, 

(-i) fc ^r fc (*) = (-ir^Ur-M, 

to account for negative angular momenta in a natural way: 

= (-l) l K-l,l- (4) 

This allows one not to use absolute value signs on /, which is quite common in the literature, but would "break" all 
the main formulae in the present work. The notational price paid is that e.g. and ^1.-1 = — ^qi do not belong 
to the same energy level: the negative-momentum state corresponding to \?xi is ^2,-1 = — ^n- 
To calculate two-body matrix elements, transform as usual to the CM frame, 

Denoting the CM basis \NL\nl) = ^> jvx(R)$ n j(r), the Coulomb interaction turns out to be 
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where Zeilberger's algorithm^ was used to reduce the customary double sum 9 to a single one, as explained in 
Appendix [X] Note that the sum has few terms in practice, since n is at most one-half the principal quantum number 
2n + I. Dropping the CM quantum numbers N,L, which introduces a notation for the relative-coordinate wave 
functions, one observes 

(n0|^|00) = (0n|^|0n), (7) 

which is quite noteworthy: the amplitude to scatter from the 2n-th level to the ground state is the same as the energy 
shift felt by the state of maximal angular momentum in the n-th level. Thus the Coulomb force has the potential to 
mix states very strongly, and this is enhanced by the cone-shape anharmonicity, as explained below. 
To get some feeling for the values involved, note that the inequality 

2n\ 4™ 

< "7= (8) 



becomes an equality in the limit n — > 00, i.e. is not an asymptotic expression. This provides the approximation 7 

(0n| ] , \0n) ~ -^=. (9) 
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The energy shifts in states of smaller angular momentum are larger, but not by much, the largest, derived in Ap- 
pendix [B] being 
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As shown in Fig. [TJ the range of values of the diagonal Coulomb matrix elements is less than a factor of two, mostly 
given by the first (square root) factor in brackets. The analytic estimates are remarkably good, all the way down to 
n = 1. 
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FIG. 1: Coulomb matrix elements (n0| |n0) (upper points) and (0n| ^7= |0n) (lower points), along with their analytic 
estimates (full lines). The intermediate case (n, n/10| |n,n/10) is also shown for n > 10 (the step at n = 20 comes from 
discrctizing n/10). All matrix elements are in units of (00 1 -±= 1 00) = y^A 



III. LABORATORY-FRAME MATRIX ELEMENTS 



Many-body wave functions in laboratory coordinates best exploit the symmetry between the particles. On the 
other hand, interaction terms in the Hamiltonian are more easily calculated in the CM frame, as above. Denoting 
the wave function in laboratory coordinates n 2 ^] = ^nih( r i)^n 2 h( r 2)i we need the transformation coefficients 

appearing in 
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They can be found following the observation^ that they must also appear in the transformation of the simplest 
monomial having the required symmetry, which is in this case 

r 2n+l e ti4> = ( x _ iy Y[ x + iy) n+l . 

One easily obtains, referring to Eq. (JSJ) , 
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where the square-root factors are needed to remove the standard normalizations of the wave functions, and Laguerre 
polynomials within them, so that the highest power of r appears with coefficient unity, for the cf s to act upon as in 
Eq. (fT2")l . The inverse of the CM transformation, Eq. (|5|), looks the same, so the same expression (|12p must hold upon 
exchanging Rf>ri and r «-» r 2 , meaning that ")" and "]" can exchange places in Eq. (TTT1) . Thus, iterating Eq. (fTTj) 
is the same as inverting it, i.e. the coefficients are idempotent: 
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Collecting the above results, one finds a computationally efficient expression for the matrix element in laboratory 
coordinates, 

[n 1 l 1 ;n 2 h\ i r \nsh; ndi] = 
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where it is assumed that An = n 3 + — n\ — n 2 > 0, without loss of generality. 



IV. THE MODEL 



The model used for shape effects is a simple modification of the single-particle spectrum. The diagonal oscillator 
energies are modified to read 

E{n,l,6) = (2n + l + l) 1+s hu, (17) 

where 6 = is of course the original oscillator spectrum, and 5 = 1 would correspond to the infinite-well spectrum, 
where the level energies increase quadratically with the principal quantum number. Thus positive values of 5 are 
referred to as "pot-shaped" dots. Negative values of d indicate level spacings decreasing with energy, and this is called 
"cone-shaped." The extreme cone shape is the hydrogen spectrum, 5 = —3 (with a negative prefactor). The idea here 
is to vary 5 in the more modest range 5 = ±0.5, with some reasonable effective Coulomb strength. 
Concerning the latter, the operative parameter is 

a = — - — , (18) 



where e is the relative dielectric constant of the parent material, and Iq = y/H/(mui) is the oscillator length scale, so 
that a measures the Coulomb repulsion at well size in units of the oscillator spacing. 

To summarize, the model considered in the present work is given by the Hamiltonian matrix 

(A\H\B) = (A\B) E ^j^+a{A\V c \B), (19) 

{n,l)&A 

where A, B are Slater determinants in configuration space and E(n, 1,5) is the modified harmonic-oscillator spectrum, 
Eq. (|17p . Calculations are performed in the subspace of minimal spin projection (zero for even number of electrons), 
which is first broken down into blocks by diagonalizing the square of the total spin operator, thus identifying states 
by spin. This achieves the same economy as if one had generated a good spin basisi^ at the beginning. 



V. TWO ELECTRONS 



Quantum dot helium, with only two electrons, is obviously not the best choice to demonstrate collective effects. 
However, it already points to the reason why cone deformations lead to enhanced collectivity. Figure [2] shows the 
evolution of the first six excited states as a function of Coulomb strength, in the usually investigated range a < 10. 

There is an obvious overall scaling effect of the anharmonicity, so to compare with the harmonic case, one should 
set energy windows to contain the same number of states, rather than the same energy range. When this is done, 
as in Fig. [2] a clear qualitative difference emerges. The cone anharmonicity induces a much greater number of level 
crossings, and they shift to lower values of the Coulomb force. 

The reason for the effect is that the higher oscillator shells open a larger phase space than the lower ones, and they 
fan out both up and down relative to the ground state. Simply enough, the states going down cross sooner with those 
coming up. Because there are many of them, the overall impact on the low-lying states is quite substantial, resulting 
in a multitude of level crossings where there have previously been only a few. 
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FIG. 2: Energy spectra of quantum dot helium, in dependence on the Coulomb strength a. Left is the recalculated result of 
Ref. H with S = in Eq. (JTTJ) , and only the first six excited states shown, the highest three degenerate. Right is the same 
with S = —0.5. All energies are relative to the ground state, which appears as a horizontal line. Successive levels in energy are 
drawn alternately by thin and thick lines. Five oscillator shells are included. Note that each "kink" on the higher curves in the 
right panel marks another level crossing, with even higher states, not shown. 
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FIG. 3: Energy spectra of four electrons for a — 3 and three values of <5, left to right: S — —0.5, 8 = 0, and S = 0.5, respectively. 
The ground state is Hund's rule triplet, denoted 3S. Note the inversion in order of the next three states in the left panel, and 
a 51 intruder, not appearing in the other two panels. The first 25 states are shown in all three panels, connected by a full line 
in order of increasing energy. Four oscillator shells were used. 




VI. FOUR ELECTRONS 



Predictably, the effect described in the previous section becomes more interesting with a larger number of electrons. 
Figure [3] shows a four-electron calculation, with 5 = —0.5, 6 = 0, and 5 = 0.5, left to right. The first thing that 
meets the eye is the loss of high orbital momentum states as 5 increases, in particular the 3H state of the oscillator 
disappears upwards for the pot anharmonicity. Conversely, a 51 intruder moves down for the cone anharmonicity, to 
become the 16th excited state, while it was 36th in the harmonic oscillator. 

Closer inspection reveals two remarkable peculiarities of the cone. First, while the 51 spin quintuplet has gone down 
in energy order, the 5S quintuplet, also present in the two other panels, has gone up, such that it is nearly degenerate 
with the 51 state of the cone. Calculation of electromagnetic transitions is beyond the scope of this work, but it is 
clear from general considerations that a 5S state lies in the decay path of a 51 state, because one can turn into the 
other by shedding orbital momentum only, without reconfiguring the spins. When they are nearly degenerate, decay 
of one into the other is likely to be significantly supressed, because of the scarcity of low-energy photons with the 
required high orbital momentum. This does not prove, without a calculation, that the 51 will have an anomalously 
long half-life, but does establish a plausible scenario for such an outcome. 

The second cone peculiarity is the inversion of order of the first three excited states, which is 1D-1S-3P for the 
harmonic oscillator, and 1D-3P-1S for the cone. The interest here is that these states have recently been the subject 
of an extensive experimental analysis^ Thus cone-anharmonic dots have a qualitative fingerprint, in principle accesible 
to experiment, which may provide some motivation to try to manufacture them, and distinguish them experimentally 
from harmonic or pot-anharmonic ones. 

The pot anharmonicity is much less interesting in comparison. Depletion of higher oscillator shell contributions 
suppresses the effects of the Coulomb force. The resulting spectra are similar to those obtained with a harmonic 
oscillator with a weaker Coulomb interaction. 



VII. DISCUSSION AND CONCLUSION 



The main result of this work is that cone anharmonicities in the single-particle level spacing of a quantum dot have 
a major impact on the collective effects, when the Coulomb force is taken into account. In the opposite limit, when 
the cone shape turns into a pot shape, the Coulomb force can be partially negated, and the spectrum recovers a form 
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similar to the unperturbed oscillator with a smaller value of the interaction. In brief, cone deformations enhance 
Coulomb collective effects, pot deformations suppress them. 

The physical mechanism leading to these effects is clear. Deviations from the single-particle spectrum, considered 
in this work, affect the higher oscillator shells more strongly. These are unoccupied in the non- interacting limit, for 
the fillings considered here. The stronger the interaction, the greater the admixture of the high shells to the ground 
state. Working against it is the energy price of involving them to optimize the occupied low level states. The benefit 
of even a slight reduction in their energy is high, because the higher shells have large degeneracies, so they make a lot 
of configuration space available for optimization. Conversely, when the component of high shells decreases, the space 
for interactions becomes more restricted, so the system behaves as if the Coulomb force were weaker than it is. 

This interpretation is confirmed when one looks at the basis-vector composition of the low cone states, such as 
the 51 state in Fig. [3] Nothing dramatic happens with respect to the analogous oscillator state, only the higher 
(in single particle ordering) components are visibly more populated, and the lower ones, visibly less. However, 
simultaneous energy and wave-function rearrangements can have an order-of-magnitude impact on decay paths and 
lifetimes, providing a mechanism to establish time-limited "quantum protectorates"— with tunable properties. 

The calculations are schematic in the sense that harmonic-oscillator wave functions do not diagonalize whichever 
potential would give the anharmonic diagonal energies, introduced by hand into the model. Nevertheless, it is not 
expected that the use of exact wave functions would affect the main outcome, because the mechanism described above 
is robust to model details. Furthermore, it is unlikely that the profile of a dot potential will ever be known with 
precision. The most that can be hoped is that some manufacturing prescription will turn out anharmonic dots with 
strongly collective charge states. If this happens, naturally one will assume that the potential well is cone-deformed. 
One might argue that the converse has already happened, and that the success of perturbation theory in describing 
real dots is partially due to the tendency of current manufacturing procedures towards pot-like anharmonicities, which 
suppress the collective effects of the Coulomb force. 

While the effect itself is robust enough, its detailed outcomes appear to be quite uncontrolled. However, in practice it 
would not be necessary to achieve highly reproducible manufacturing to observe the effects discussed above. (Quantum 
dot molecular arrays are currently also produced with a considerable range of properties^) All that would be required 
to first order is that an acceptable percentage of manufactured dots have a tendency to "get stuck" in their de-excitation 
pathway. In any case, it is far-fetched at present to speculate on the possible ramifications of the simple proof of 
principle, given here. It is hoped, on the contrary, that its simplicity and universality will motivate someone to try 
to obtain strongly collective few-electron states by manufacturing anharmonic quantum dots with cone-like confining 
potentials. 



Appendix A: Coulomb interaction in relative coordinates 

The usual expression^ for the Coulomb matrix element may be rewritten 



where 
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Zeilberger's algorithm 5 gives^ 
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which could in principle be checked directly. The solution is immediate, 
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Inserting it above and rearranging terms into binomial coefficients gives Eq. In checking the latter against the 
literature, one should take into account the caution below Eq. Q. 



Appendix B: Approximation for the zero-angular momentum matrix element 

Using the approximation Eq. flS)) for the case I — in Eq. (J6j> , one easily finds 
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where the end terms in the sum are treated specially to preclude a nonsensical result. The sum is approximated by 
an integral, 

I - "" 1 1 , 21n(n-l) 
/ —, ^dx = , 

J 1 x(n — x) n 

and finally ln(n — 1) is replaced by Inn by inspection, because it gives a much better approximation. 
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